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ABSTRACT 


The  purpose  of  this  study  was  to  analyze  and  compare  papers 
by  Lee  (ref.  7),  Thompson  (refs.  9  and  10)  and  Ingold  (ref.  6) 
which  give  conflicting  results  concerning  the  charged  particle 
profiles  in  an  oxygen  discharge  tube.  Lee  predicts  proportional 
negative  ion  and  electron  profiles  whereas  Thompson  and  Ingold 
predict  non-proportional  profiles. 

The  analytic  developments  were  critically  reviewed  and 
numerical  solutions  were  developed  for  each  approach.  It  was 
found  that  Thompson's  results  could  not  be  obtained  without 
introducing  the  additional  assumption: 

V«e 

Using  this  relation  and  the  assumption  that  a* i  is 

indicated  by  Thompson's  experimental  results)  a  new  analytic 
relation  was  developed  for  the  negative  ion  profile  which  agreed 
well  with  Thompson's  experimental  results. 

Ingold's  development,  on  the  other  hand,  seemed  to  be 
lacking  a  firm  mathematical  basis.  One  of  Ingold's  fundamental 
relations  used  throughout  his  development  could  not  be 
established  and,  furthermore,  it  was  shown  that  a  completely 
different  solution  could  be  developed  without  introducing 
further  constraints.  Also,  the  numerical  solution  developed 
for  Ingold  did  not  give  profiles  similar  to  Ingold's. 


Lee's  profiles  reduce  to  Thompson's  when  associative 
detachment  is  ignored  and  a  common  set  of  parameters  and 
coordinate  system  are  used.  When  associative  detachment  is 
included,  Lee  and  Thompson  do  not  give  similar  profiles.  Thus, 
the  relation 


may  only  be  valid  in  limiting  cases  (e.g.  when  associative 
detachment  can  be  ignored.)  Results  from  an  investigation  of 
power  series  solutions  indicated  that  they  could  probably  be 
used  to  obtain  Lee's  profiles,  thus  negating  the  need  to 
numerically  integrate  simultaneous  first  order  equations,  as  is 
done  by  Lee. 


ix 


I.  INTRODUCTION 


The  low  pressure  gas  discharge  with  two  charged  species, 
positive  ions  and  electrons,  was  modeled  by  Schottky  in  1924. 
Predictions  of  his  theory  have  agreed  well  with  experimental 
results.  Recently,  though,  attention  has  been  focused  on  the 
behavior  of  electronegative  gas  discharges  because  of  their 
applicability  to  semiconductor  processing,  lasers  and  ion 
sources.  To  model  electronegative  gases,  Schottky  theory  must 
be  extended  to  include  negative  ions.  Several  attempts  have 
been  made  to  incorporate  negative  ions  but  results  have  not  been 
in  agreement.  Lee  claims  that  the  charged  particle  profiles  are 
proportional  to  one  another  but  Thompson  and  Ingold  claim  the 
profiles  are  not  proportional. 

Problem 

This  study  will  critically  analyze  current  papers  on  three 
component  diffusion,  outline  assumptions,  and  develop  numerical 
solutions  in  order  to  compare  and  contrast  the  various 
approaches  and  possibly  draw  conclusions  concerning  their 
validity. 

The  study  will  be  limited  to  papers  by  Thompson  (Ref  9  and 
10) ,  Lee  (Ref  7) ,  and  Ingold  (Ref  6) . 

Assumptions 

All  papers  studied  begin  with  the  diffusion  equations 

(1-la)  II  "  E 


1 


r.=  -D.Vn.  E 

rt° 


(1-lb) 

(1-lc) 


where 


is  the  particle  flux  density  and  D  is  the  diffusion 


coefficient  which  can  be  defined  through  the  Einstein  relation 


(1-2) 


A  is  the  mobility  defined  as  (4:186) 

ju  =  -3i- 


(1-3) 


where  is  the  collisional  frequency  between  the  species  and 

the  neutral  particles  in  the  plasma,  n  represents  the  species 
number  density  and  E  is  the  electric  field.  Thus,  it  is  assumed 
that  the  diffusion  equations  are  applicable.  The  assumptions 
necessary  in  obtaining  the  diffusion  equations  from  first 
principles  are:  no  external  forces  acting  on  the  charged 
species  other  than  an  electric  field,  neutral  particles  have  no 
mean  velocity,  pressure  is  hydrostatic,  and  the  drift  velocity 
of  the  charged  species  is  assumed  to  be  negligible  compared  to 
the  mean  thermal  velocity. 


Research  began  with  a  thorough  analysis  and  critique  of  each 
paper  studied.  This  was  necessary  before  any  further  work  could 


be  initiated  and  offered  some  insight  as  to  how  much  further 
work  each  paper  warranted.  Next,  numerical  solutions  yielding 
particle  profiles  were  developed  for  each  approach.  This 
facilitated  comparison  and  contrast  especially  when  changes  and 
variations  were  made  to  each  theory.  Next,  an  unsuccessful 
attempt  was  made  to  directly  integrate  the  three  diffusion 
equations,  three  continuity  equations  and  Poisson's  equation 
(seven  equations  and  seven  unknowns) .  Finally,  a  more  detailed 
analysis  was  made  of  the  power  series  expansions  (both  on  axis 
and  at  the  wall)  used  by  Lee  to  investigate  the  possiblities  of 
using  power  series  expansions  to  obtain  profiles  over  the  entire 
tube  radius,  thus  eliminating  any  need  for  direct  integration  of 
equations. 

This  presentation  will  begin  with  the  assumptions  made  by 
each  paper,  followed  by  discussions  of:  the  derivation  of  the 
equations  found  in  each  paper  (many  of  these  derivations  can  be 
found  in  the  appendices) ,  the  developments  of  the  numerical 
solutions,  the  behavior  of  the  equations  in  the  two  component 
limit  and  the  sensitivity  of  each  solution  to  a  variation  of 
parameters.  The  presentation  concludes  with  a  comparison  of 
results  and  a  short  discussion  of  boundary  conditions. 

Backjcound 

As  important  backround  information,  we  will  now  briefly 
redevelop  the  two-component  Schottky  solution  (4:187). 


We  assume 


where  n  stands  for  charged  particle  density  and  J»  is  particle 
current  density  or  flux.  1-1  can  then  be  written  as 


f  a-kai 


where  the  gradients  are  in  Cartesian  coordinates.  Eliminating  E 
and  solving  for  £  ,  we  have 


(1-6) 


where 

in  -  +A-  P+  _  + 


(1-7) 


This  last  relation  was  obtained  by  employing  the  Einstein 
relation 

Oj  .  k£  =  § 

A*'  t  ~  1  (1-8' 


Using  the  continuity  equation 


II.  ASSUMPTIONS  AND  ANALYTIC  DEVELOPMENTS 


Besides  the  assumptions  implicit  in  the  diffusion  equations^ 
certain  other  assumptions  are  made.  Three  of  these  are  common 
to  all  three  papers  analyzed.  The  first  of  these  is 
quasineutrality  i.e. 


ft-  +  (2-1) 

The  second  is  that  the  net  flow  of  charge  is  zero  everywhere 
i.e. 

(2-2) 

The  third  is  that  both  the  negative  ion  and  positve  ion 
temperatures  are  equal  to  the  backround  gas  temperature. 

Lee  further  assumes  that  there  are  three  production  and  loss 
terms: 

ionization*  *0  (T«)n,cnf 3  sec”1  (2-3a) 

dissociative  attachment*  Afc)  cm”  3  sec-1  (2— 3b) 

associative  detachment*  (pCMt'k  cm  sec  (2-3c) 

Lee  also  assumes  that  the  charged  particle  densities  are  zero  at 
the  outer  boundary  and,  for  his  analytic  solution,  he  requires 
the  following  proportionality  relation  to  be  true: 


Ingold  uses  ionization  and  dissociative  attachment  as 
production  and  loss  terms  but  ignores  associative  detachment 


Also,  like  Lee,  Ingold  assumes  the  particle  densities  go  to  zero 
at  the  outer  boundary. 

Thompson  uses  the  same  production  and  loss  terms  as  Ingold. 
Furthermore,  although  Thompson  doesn't  specifically  state  the 
assumption  of  a  proportionality  relation,  it  was  determined  that 
the  relation 


Vn- 

?n«. 


- 


(2-5) 


is  needed  to  obtain  Thompson's  results.  (See  appendix  A) 
Discussion  of  Derivation  of  Equations 

Like  Edgley  and  Von  Engel  (2) ,  Lee  initially  develops  a 
system  of  seven  equations  and  seven  unknowns.  These  seven 
equations  are  the  three  diffusion  equations,  the  three 
continuity  equations  (see  eqs.  B-l  and  B-2)  and  Poisson's 
equation 

V-E  = (2_( 


These  are  reduced  to  a  system  of  four  equations  and  four 
unknowns  by  imposing  2-1  and  2-2  and  eliminating  the  electric 
field  from  the  set  of  equations.  In  general,  Lee's  step  by  step 
development  is  detailed  enough  such  that  analysis  of  his  work 
consisted  mainly  of  checking  for  algebraic  mistakes  —  none  of 
which  were  found  in  his  analytic  development.  Thus,  no 


appendix  was  necessary  for  tracing  the  derivation  of  Lee's 
equations. 


Ingold,  like  Lee,  eliminates  the  electric  field  from  his 
equations.  His  development,  however,  seems  to  be  seriously 
flawed  in  that  the  relation  (eq.  B-26) 

♦•k'y  )=  y +  )  (2-7) 

which  is  essential  to  his  derivation,  could  not  be  established 
(See  appendix  B.)  Rather,  the  relation  (eq.  B-30) 

(2-8) 


where  a,b,c  and  d  are  constants,  seems  to  be  the  correct 
equation. 

It  is  possible  to  re-develop  Ingold's  equations  using  2-8. 
The  derivation  is  as  follows. 

/ 

Dividing  2-8  through  by  X  we  have 

(2-9) 

As  the  spatial  coordinate,  z,  goes  to  zero,  the  ratio 
becomes  indeterminate  since  both  derivatives  go  to  zero.  Usin 
L'hopital's  rule  we  obtain 


where  the  subscript  indicates  the  variable  is  evaluated  at  z**0. 
Thus  2-9  becomes 


(2-11) 


Multiplying  through  by  and  rearranging,  we  have 


(2-12) 


u 

Dividing  through  by  yo  ,  we  obtain  a  quadratic  equation  with 
the  solution 


\ 


c-L± 

2«. 


(2-13) 


This  relation  replaces  B-29  and  must  be  used  in  it's  place  when 
substituted  into  B-24.  Thus  B-24  becomes 


e*x. 


«  v  a  f  lfe‘.  V 

■  ^  W  ♦  i  ♦  1 


■] 


=  e>y. 


A. 

A' 


*r;j  ? 


•*  A* 


/*■ 


n«o 

n.o 


(2-14) 


which  simplifies  to 


(2-15) 


which  differs  considerably  from  B-31 


The  next  change  necessary  is  to  replace  B-44  with 


(2-16) 


With  this  change,  we  can  proceed  exactly  as  Ingold  does  to 
obtain  the  parametric  solutions 


y*C'uS*&«4 


(2-17a) 

(2-17b) 


as  opposed  to  B-50  and  B-52.  The  difference  between  A  and  A' 
and  C  and  C'  is  that  A  and  C  can  be  determined  by  the  boundary 
conditions 


Xo 


(2-18) 


whereas  A'  and  C'  cannot  be  determined  by  the  parallel  boundary 
condition 


XJ  -  0  S/J  -  0  (2-19) 

Therefore,  another  boundary  condition  would  be  needed  for 
further  development  of  this  solution.  If  B-54  is  differentiated 


then  2-17a  and  b  can  be  substituted  in  to  give  an  equation 
identical  to  B-56  except  that  A  and  C  are  repaced  with  A'  and 
C' . 


Even  this  solution  seems  somewhat  arbitrary 
demonstrated. 

It  was  pointed  out  in  appendix  B  that  eq.  B 
an  arbitrary  assumption  fixing  the  value  of  f. 
fix  the  value  of  f  such  that 

k-lf-0  or  t--  Va 

Then  B-49  becomes 


du  _  V\  +  |0  M 
dix'  ***' 

where  x  is  replaced  by  x*  to  conform  with  2-16. 
constants  take  on  new  values  such  that 

)rr\-  at  b'f 
H  -  CtoK  j 


as  will  now  be 

47  seemed  to  be 
Suppose ,  we  now 

(2-21) 


(2-22) 

Also,  the 

(2-23) 


Rearranging  2-22,  we  have 


Making  the  following  definitions 


Y=R\f£  ^ 

and  substituting  2-30  and  2-27  into  2-20  we  obtain,  after 
simplification. 


=o 


(2-32) 


The  solution  to  this  second  order  differential  equation  can  be 
seen  by  inspection  to  be  some  combination  of  cosines  and  sines. 
This  differs  considerably  from  B-55  which  has  no  obvious 
analytic  solution. 

Thompson,  like  Lee  and  Ingold,  eliminates  E  from  his 
equations.  Furhtermore,  an  unusual  proportionality  relation 


El- 

Vfle 


(2-33) 


is  needed  to  obtain  Thompson's  ambipolar  diffusion  coefficients 
(see  appendix  A) .  The  fact  that  this  relation  is  required  can 
be  further  demonstrated  by  substituting  A-2  into  the  assumption 
of  zero  total  charge  flux 


-D>*7a»  =  -Dlawl~D*  Vn*. 


(2-35) 


Since  7 »4*VA-+7h«  ,  we  can  simplify  2-35  to 

Vn.  ,  -  0 

?rte '  D* - b/ 


(2-36) 


i 


From  A-3,  we  can  evaluate  the  right  side  of  2-36  to  give 


1  +  */»-/*»«  _  \ 

\  »Y«<  » 

-*■  At  -  I 
r  <*•  L  **<Y 


-  HiL&iP.wr 

»  a!  "• 


(2-38) 


Therefore 


VW; -  _  ,*y 

'Jd.  (2-38) 


This  relation  was  required  to  obtain  Thompson’s  expressions 
for  D£  ,  D*  and  D*  however  D«  could  also  be  obtained  from  the 
assumption  that  . 

Using  Cartesian  coordinates,  as  Thompson  does,  we  can  write 
2-38  as 


14 


(2-39) 


*• 


1  ^L“  ,yi  4^ 
h-  X£  cTS 


Integrating,  we  have 


j2/v\  "*  H-o  ”V f " Xv>. ^ ^  "t  C. 


(2-40) 


Taking  the  exponential  of  both  sides  we  obtain 


^  -cf^Y 

ru  o  ^WoJ 


(2-41) 


where  C=1  since 


l*C(0  *t  *~-0 


(2-42) 


An  interesting  relation  can  be  developed  if  we  assume  j^a ^  | 
as  is  indicated  by  Thompson's  experimental  results  throughout 
most  of  the  discharge  tube.  Taking  the  power  series  expansions 
for  ln(l+x) ,  we  have 


(2-43) 


If  we  define  x  as 


-  Hi  -  ) 
‘  n«0  * 


(2-44) 


15 


t.i 


£| 

—■-i 


then  x<<1  and  we  can  write 


JU(l*  te-O-JU^.  - 


(2-45) 


So  2-40  becomes 


or 

n-  -  0  (2_46) 

where  Vt  -electron  temperature  and  V-  -negative  ion 
temperature.  This  resembles  the  Boltzmann  distribution  given  by 
Thompson  (10:820)  to  model  the  negative  ion  concentration  as 
given  by  his  experimental  results.  The  exact  form  given  by 
Thompson  is 

]  (2-47) 

where  V  is  the  spatial  potential  relative  to  the  discharge  axis 
(z«0) .  On  page  819  of  ref.  10,  Thompson  plots  v/V,  according 
to  experimental  measurements. 

From  Thompson's  graph  of  (10:820)  one  can  see, 

without  too  much  imagination,  that  a  graph  of  (  -  1) 

would  look  very  similar  to  a  graph  of  V.  It  is  difficult  to  say 


how  closely  they  match  since  Thompson  doesn't  give  any  specific 
numbers  and  his  graphs  obviously  weren't  drawn  for  the  purpose 
of  reading  off  specific  data  points.  Nonetheless,  it  is 
interesting  to  note  that  an  analytic  expression  can  be  obtained 
from  the  assumptions 


y  and  Hj; 


which  seems  to  agree  fairly  well  with  experimental  results  given 
by  Thompson . 


Because  Lee  uses  cylindrical  coordinates,  he  is  forced  to 
deal  with  singularities  at  the  wall  and  on  axis.  The 
singularity  on  axis  is  dealt  with  by  using  the  power  series 
expansions  given  by  06.  The  coefficients  are  found  through  the 
recursion  relations  Oil,  014,  016,  and  023  and  then  a  small 
value  for  z  is  chosen  so  initial  values  for  y,  ,  yr  ,  y  a  and 
can  be  obtained  for  use  in  the  integration  of  equations  01a,b,c 
and  d.  Fig.  1  depicts  the  normalized  charged  particle  profiles, 
given  by  the  code  outlined  in  appendix  D,  as  functions  of  the 
normalized  radial  coordinate,  z. 

The  electric  field,  as  evaluated  in  the  code,  is  shown  in 
fig.  2.  The  field  is  seen  to  rise  sharply  as  one  leaves  the 
axis,  pass  through  an  inflection  point  at  about  z*0.5  and  rise 
asymptotically  as  the  outer  wall  is  approached.  This  latter 
behavior  is  to  be  expected  since  we  have  required  the  charged 
particle  densities  to  go  to  zero  at  the  wall. 

The  singularity  at  the  outer  wall  is  dealt  with  by  the 
formal  approximations  given  by  C-24.  It  is  shown  in  appendix  C 
that  C-24  reduces  to  C-25  when  terms  of  order  and  larger 

are  neglected.  The  integration  routine  proceeds  to  a  point,  z, 
a  small  distance  ,  ,  away  from  the  outer  wall.  The  values 

of  y,  and  y->.  at  this  point  are  used  to  find  and  L>(  in 

C-25.  Given  and  b4  ,  we  can  find  C0  and  as 


FIG.  1  CHARGED  PARTICLE  PROFILES  (LEE) 


defined  in  C-31  and  C-33.  Given  Ca  end  d0  ,  we  are  now  able 
to  evaluate  y a  and  ^4  as  given  by  C-25  (at  1-6,.  )  and 
compare  these  values  with  those  obtained  from  the  integration 
routine  at  the  point  1*6,.  .  If  the  two  sets  of  values  differ, 
then  the  adjustable  parameters  —  pressure,  electron 
temperature,  and  h  —  must  be  varied.  The  profile  in  fig.  1 
was  obtained  by  holding  the  pressure  constant  at  0.4  torr  and 
varying  h  until  a  minimum  was  found  in  the  differences  between 
the  values  of  ^3  and  obtained  from  the  integration 

routine  and  the  Taylor  series  expansion  near  the  wall.  The 
electron  temperature  was  then  varied  and  the  process  was 
repeated.  This  was  done  a  number  of  times  to  find  the  best 
electron  temperature  for  the  given  pressure.  Further  discussion 
of  the  sensitivity  of  the  profiles  to  a  variation  of  the 
parameters  will  be  made  at  the  end  of  the  chapter. 

Finally,  although  Lee  intended  his  power  series 
approximations  on  axis  to  be  used  only  in  the  near  vicinity  of 
the  axis,  it  was  found  that  the  series  gave  profiles  which 
agreed  remarkably  well  with  those  obtained  from  the  numerical 
integration  of  the  equations  up  to  values  of  z  greater  than 
0.9.  Figs.  3a  and  3b  show  profiles  obtained  from  both  the  power 
series  approximations  and  the  integration  routine.  3a  is  a  plot 
of  normalized  electron  densities  versus  normalized  spatial 
coordinate,  z,  while  3b  depicts  normalized  negative  ion 
densities.  The  first  100  terms  in  the  power  series  expansions 
(eqs.  C-6)  were  used  for  the  series  approximations.  As  is 


evident  from  the  graphs/  excellent  agreement  with  the 
integration  routine  was  obtained.  Figs.  4a  and  4b  demonstrate 
the  convergence  behavior  of  the  power  series  approximations  at 
z*0.9  and  z=0.95.  The  magnitude  of  the  term  in  the  power 
series  was  plotted  as  a  function  of  k.  In  both  cases/  the 
magnitudes  decreased  rapidly  for  k  greater  than  approximately  90 
indicating  convergence  had  been  reached.  The  series  does  not 
converge  at  z«l;  thus,  the  power  series  approximation  begins  to 
diverge  between  z«.95  and  z«l. 

A  fourth  order  Runge-Kutta  method  (5:95)  is  used  for  the 
numerical  integration  (solution  of  a  second  order  differential 
equation).  Ingold's  parametric  solution  as  expressed  in 
equation  B-56  cannot  be  solved  using  a  Runge-Kutta  integration 
technique  since  the  technique  can  only  handle  first-degree 
differential  equations  (i.e.  the  derivatives  are  not  squared, 
cubed,  etc.)  This  problem  can  be  overcome  by  recognizing  the 
fact  that  B-56  can  be  written  as 

Bu)  -  0  (3-1) 


where 


W  =■  (A  ♦  ^(C)  US+(&+  D)  <4 


A  new  problem  is  encountered  in  that  3-1  now  contains  two 
variables.  This  problem  is  dealt  with  by  using  an  iterative 


subroutine  that  finds  the  appropriate  value  of  u  for  a  given 
value  of  w  using  equation  3-2.  Thus,  the  subroutine  is  called 
whenever  a  value  of  x  is  needed  where  the  last  value  obtained 
for  w  is  used  to  find  u. 

Mobilities,  ionization  rates  and  dissociative  attachment 
rates  are  taken  from  Thompson  since  Ingold  doesn't  supply  any  of 
these  values.  The  numerically  generated  profiles  obtained  from 
the  code  listed  in  appendix  F  did  not  resemble  Ingold's.  Rather 
than  obtaining  non-proportional  profiles  as  indicated  by  Ingold, 
the  profiles  were  seen  to  be  proportional  over  a  wide  range  of 
values  for  h. 

Thompson's  only  analytic  expressions  are  those  he  obtains 
for  his  ambipolar  diffusion  coefficients.  However,  it  is  shown 
in  appendix  A  that  the  proportionality  relation 

—  =  Y°< 

V». 

is  necessary  to  obtain  Thompson's  diffusion  coefficients.  From 
2-41  we  have 


Equation  B-9  is  a  general  relation  for  which  no  assumptions 
specific  to  Ingold  were  made,  thus,  we  will  use  it  here  and 
employ  3-3  to  obtain 


(3-4) 


[yCy-O^'W  YxtV>  k'x  -  0 


As  with  Ingold,  this  can  be  written  as 


W7/  +  x  =  0 


(3-5) 


where 


W  = 


e.+Sj  * 


+  X 


(3-6) 


The  coding  of  these  equations  for  numerical  integration 
proceeded  exactly  as  Ingold's.  An  iterative  subroutine  was 
called  when  necessary  to  find  the  appropriate  value  of  x  for  a 
given  w  using  3-6. 

Again  a  fourth  order  Runge-Kutta  (5:95)  method  was  used  for 

the  numerical  integration.  Thompson  uses  a  Boltzmann 

distribution  to  model  the  negative  ion  density  and  concludes 

that  a  negative  ion  temperature  of  0.15  ev  gives  the  best 

agreement  between  theoretical  results  and  experimental  data; 

thus,  0$  was  taken  to  be  0.15  ev  in  the  numerical  code. 

Thompson  also  gives  V  =16  which  fixes  ©  e,  at  about  2.4 

2  —1  —1 

ev.  is  given  as  2.25  cm  V  sec  (9:517)  and 

Thompson  states  that  *0.0043  and  =0.0022  for 

2  -1 

oxygen.  This  fixes  yt»«.  at  1022  and  at  4.39  cm  V 

sec-*.  Using  these  values  for  the  mobilities  and  ion 


temperatures/  it  was  found  that,  to  obtain  profiles  similar  to 


those  given  by  Thompson,  ionization  and  dissociative  attachment 

-1 


rates  on  the  order  of  10  sec  "*■  were  needed.  The  ionization 
rate  is  the  same  order  of  magnitude  as  that  predicted  by  the  two 
component  Schottky  limit  where 


fa?  „  3L 
V  -o  “  z 


(3-7) 


The  above  relation  yields  a  value  of  *0  =14.1  sec-^  for  the 

mobilities  and  ion  temperatures  given  by  Thompson.  If 
ionization  and  dissociative  attachment  rates  similar  to  Lee  are 
used  ( 1 0 6  and  105  sec  1) ,  then  mobilities  on  the  same 
order  of  magnitude  as  Lee  (10^  for  electrons,  10^  for  ions) 
are  needed  to  again  obtain  profiles  similar  to  Thompson's. 

Exact  values  for  "0  and  X  are  difficult  to  obtain  since 
Thompson  makes  no  attempt  to  establish  any  sort  of  boundary 
conditions  and  none  are  apparent  upon  inspection  of  his 
profiles.  Figs.  5  and  6  show  several  profiles  for  various 
values  of  "0  and  X 
Two  Component  Limit 

A  requirement  that  all  solutions  to  the  three  component 
problem  should  satisfy  is  that  they  give  results  similar  to  the 
analytic  two  component  solutions  (where  only  electrons  and 
positive  ions  are  considered)  in  the  two  component  limit.  This 
limit  can  be  obtained  in  Lee's  code  (appendix  D) ,  if  two 
variables  ,  h=  and  alpha,  are  reduced  to  small  values. 


NORMALIZED  DENSITY 


NORMALIZED  DENSITY 


•  M 


32 


NORMALIZED  DENSITY 


We  can  set  h  to  zero  but  not  alpha  since  division  by  zero  would 
occur  in  a  number  of  places.  This  problem  can  be  overcome  by 


redefining  Vi  and  as 

-  iiilii  v  =  H£l> 

rather  than  the  definitions  given  by  C-3.  This  was  done  and  the 
appropriate  changes  were  made  in  the  code  and  it  was  observed 
that  the  profiles  did  not  differ  significantly  from  those 
obtained  by  using  Lee's  definitions  and  a  small  value  of  alpha 
instead  of  zero.  (Values  from  0.1  to  0.001  sec-1  gave 
essentially  the  same  profiles.) 

Since  Lee  uses  cylindrical  coordinates,  the  electron  profile 
would  be  expected  to  look  like  a  zeroth  order  Bessel  function. 

We  have  shown  that,  in  cartesian  coordinates,  the  relation 


must  be  obeyed  if  the  charged  particle 
zero  at  the  wall.  Similarly,  from  our 
obtained 

r 

and 


densities  are  to  go  to 
derivation  of  3-9,  we  had 


(3-10) 


y.f1  -  -0  v\ 


(3-11) 


where,  in  cylindrical  coordinates,  we  obtain 


(3-12) 


or 


Z  7i> 


(3-13) 


Thus,  if 


D* 


(3-14) 


the  solution  of  3-13  is  a  zeroth  order  Bessel  function.  3-14 
now  specifies  "0  given  the  mobilities  and  ion  and  electron 
temperatures  since 


-0  -  D«,= 


j.  j_ 

/*K 


(3-15) 


Besides  including  this  condition  in  the  code,  one  other  change 
must  be  made.  The  tube  radius  was  assumed  to  be  1  cm  so  that 
the  R  term  in  3-4  could  be  neglected.  We  are  now  forced  to 
include  it  and  set  it  equal  to  2.405  since  this  is  where  the 
first  zero  of  the  zeroth  order  Bessel  function  occurs.  With 
these  changes,  it  was  found  that  the  electron  profile  given  by 
the  code  agreed  reasonably  well  with  a  Bessel  function.  The 
code's  profile  did  diverge  slightly  from  the  Bessel  solution. 


increasing  to  a  maximum  difference  of  about  0.05  at  about  x»0.9 
and  then  converging  back  to  the  Bessel  solution  since  both 
profiles  go  to  zero  at  the  wall. 

To  see  how  Ingold's  equations  reduce  in  the  two  component 
limit,  we  begin  by  dividing  a,  b,  c  and  d  by  I1).*  and  redefining 
these  terms  as  a,  b,  c  and  d  and  employing  Ingold's  relation 
B-31  to  obtain 


•  fj  ( l 


CkMt 


(3-16) 


where,  again,  ft.s  and  «<,  have  been  set  to  zero.  From  B-47  we 
have 


r  ^  Ao 

f  -cTI,  (3-17> 

and 

hn  -r  =  j 

Y>  -  do  -  o  ; 

Then,  from  B-51  and  B-53  we  have 

A  ~  sA  o  j  B  -  £rn  ~  o 

c  -A4w.3t.sCo :  t>~  i  ^  f  - 1 


Y\  -  Catdsf  - 

c  r  *2  =  £2  -  C 
r*\  cTo  "  * 


(3-18) 


So  B-56  becomes 


A*  s6u^VV^*(s.-i)u^VVk'A>us*^  0 


or 

^(A.us‘J  +k"A.us-  =  o 


The  solution  of  interest  to  this  equation  is 

Ku°  *  c&s(Jc 


(3-20) 


(3-21) 


and  from  B-50,  we  see 


x  -  A,  cosQc^) 


(3-22) 


From  B-55,  we  have 


Lx  - 

^->0  *  " 


3l  +  £- 

A*  T  Af 


s 


-J 

p* 


(3-23) 


which  is  the  expected  result.  Thus,  Ingold's  analytic  equations 
reduce  properly  in  the  two  component  limit. 

Likewise,  the  code  in  appendix  F  gives  the  proper  cosine 
curve  for  the  electron  density  when  H.#  and  «<  go  to  zero  and 


the  relation 


is  obeyed,  even  though  the  code  does  not  give  profiles  similar 
to  Ingold's  in  the  general,  three  component  case. 


For  Thompson,  we  see  from  3-4  that  as  h  goes  to  zero,  we 
obtain 


(3-24) 


The  appropriate  solution  of  which  is 


Y  -  coslk*) 


(3-25) 


where  the  k  expressed  here  is  the  same  as  Ingold's  k.  Thus,  3-4 
also  has  the  proper  behavior  in  the  two  component  limit. 

We  now  turn  our  attention  to  Thompson's  analytic 
expressions.  As  x  goes  to  zero,  equation  A-3a  becomes 

it  r 


D1'-  h  r- 

“  ***■  l  t 


s  ■ 


tty 
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(3-26) 
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Ar  At 


where  we  have  employed  the  Einstein  relation.  We  recognize  3-26 
as  the  ambipolar  diffusion  coefficient  for  a  two  component 
plasma.  Also,  for  A-3c 


'  A«  T  A* 

Furthermore,  one  would  expect,  in  the  limit  as  goes  to 

infinity,  that  both  A-3a  and  A-3b  would  reduce  to  D+  if 
From  A-3a,  we  have 


(3-28) 


and  from  A-3b,  we  have 


fc*  rj)  1  hr  _ ik 

<<-joo  ^  A« 


oc  y 
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A*  At 
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(3-29) 


Thompson  also  gives  an  analytic  expression  for  the  electric 
field  as  follows  (10:20): 

E  0*0  =  1  ~  £»/&♦?  \ yV 

E(g)  ~  v  (3"30) 


where  E{0)  is  the  electric  field  for  =  0  .  Thus,  if  we  look 
at  3-30  in  the  limit  as  •<  goes  to  zero,  we  would  expect  the 
right  side  to  go  to  one.  Doing  this,  we  obtain 


A~,  IW. 

<-90  g&)  ,  '  V  / 


v  y 


(3-31) 


Thus,  Thompson's  expression  for  the  electric  field  does  not 


reduce  properly  in  the  two  component  limit.  Failing  to  meet 
this  basic  requirement,  we  shall  conclude  that  Thompson's 
electric  field  expression  is  incorrect.  An  attempt  was  made  to 
derive  a  correct  analytic  expression,  but  none  could  be  found 
that  did  not  include  a  term  similar  to  • 

Sensitivity  to  Parameters 


Lee's  profiles  were  seen  to  be  very  sensitive  to  h,  the  on 
axis  ratio  of  negative  ions  to  electrons.  Many  difficulties 
were  encountered  in  trying  to  find  a  value  of  h  which  gives  the 
best  fit  to  the  boundary  conditions.  A  measure  of  the  fit  to 
the  boundary  conditions  at  the  wall  is  given  by  the  terms  D1F1 
and  D1F2  as  found  in  the  code  outlined  in  appendix  D.  D1F1 
represents  the  difference  between  y3  ,  as  found  by  the 
integration  routine  and  yj  ,  as  given  by  025.  D1F2 
represents  the  difference  between  y^  •  as  found  by  the 
integration  routine  and  y^  ,  as  given  by  C-25.  The  fit  is 
made  at  a  point  a  small  distance,  ,  away  from  the  wall. 

Thus,  h  must  be  chosen  so  as  to  minimize  D1F1  and  D1F2. 

Because  of  the  poor  convergence,  Lee's  graphs  of  pressure 
versus  electron  temperature  and  h  versus  pressure  (L:4702)  were 
used  to  obtain  initial  values.  It  was  found  that  if  h  was  off 
by  more  than  several  tenths,  yx  would  begin  to  rise  before  the 
wall  could  be  reached.  Once  a  range  for  h  was  found  where 
was  always  negative,  the  values  where  DlFl  and  D1F2  reached 
minimum  values,  or  changed  sign,  could  be  investigated  with  a 
smaller  &  k  .  To  complicate  the  problem,  though,  convergence 


was  not  steady  and  several  apparent  minima  often  occurred,  even 
with  spacings  as  small  as  0.01.  Upon  closer  investigation  of 
some  minima,  using  a  smaller  ,  unusual  behavior  such 

as  going  positive  would  occur.  Well  behaved  convergence 

could  only  be  achieved  when  the  initial  guess  wasn't  more  than 
about  0.1  away  from  the  final  value. 

For  P*0.4  torr,  it  was  found  that  *2.1  eV  and 

h*7.2735.  This  gave  minimum  values  of  D1F1  and  D1F2  of  -0.879 
and  0.587  respectively.  Increasing  h  by  0.01  causes  DlFl  and 
D1F2  to  rise  to  -2.253  and  1.288  and  increasing  h  by  0.1  causes 
D1F1  and  D1F2  to  jump  to  -40.508  and  23.206. 

Increasing  by  0.01  eV  resulted  in  an  increase  in  h  to 

7.4743  with  DlFl  and  D1F2  remaining  at  about  the  same  values  as 
when  Ge  *2.1  eV.  Increasing  ©e  by  0.1  eV  to  2.2  led  to 
minimums  of  DlFl  and  D1F2  on  the  order  of  LO *“ 

Finally,  pressure  was  varied  and  it  was  found  that  smaller 
values  of  DlFl  and  D1F2  could  be  obtained  for  P*0.47  torr 
when  Ge  -2.1  eV.  It  was  found  that  D1F1=0.0135  and 
D1F2— 0.0083  for  h*7.4199.  It  is  possible  that  DlFl  and  D1F2 
could  be  reduced  to  even  smaller  values  given  the  right 
combination  of  pressure,  electron  temperature  and  h,  but,  given 
the  difficulty  of  trying  to  simultaneously  minimize  three 
parameters  and  the  fact  that  h  does  not  converge  steadily,  time 
limitations  prevented  further  investigation. 

Figures  5  and  6  demonstrate  the  sensitivity  of  Thompson's 
profiles  to  changes  in  i)  and  \  .  One  might  expect  that 


changes  in  ©e  would  lead  to  steeper  negative  ion  profiles 
since  -y  x,  ,  but,  it  was  observed  that  an  increase  in  ©e 
resulted  in  flatter  electron  profiles  and  relatively  no  change 
in  the  negative  ion  profiles.  Even  with  ©c  increased  to  10 
times  the  value  used  by  Thompson  and  the  code  in  appendix  E,  it 
was  seen  that  the  negative  ion  profile  changed  by  less  than  0.01 
where,  for  *0  -22  and  X  =11,  the  electron  profile  was 

almost  flat,  reaching  a  minimum  of  about  0.98  at  the  wall. 
Decreasing  resulted  in  both  steeper  electron  and  negative 

ion  profiles.  (For  ©^  =0.05  eV,  x  went  to  zero  at  0.88  and  y 

went  to  zero  at  0.61.) 

The  sensitivity  to  changes  in  parameters  in  the  code  for 
Ingold  was  not  investigated  since  it  would  have  been  pointless 
considering  all  the  problems  previously  indicated  in  its 
developement. 


IV.  RESULTS- AND . CONCLUSIONS 


Comparison  of  Results 

As  was  mentioned  earlier,  if  Thompson's  data  for  mobilities, 
electron  and  ion  temperatures  and  ratio  of  negative  ion  density 
to  electron  density  on  axis  are  used  in  the  code  for  Thompson, 
it  is  found  that  ionization  and  dissociative  attachment  rates  on 
the  order  of  10  are  needed  to  give  profiles  that  are  similar  to 
those  given  by  Thompson's  experiemental  data  (10:820).  Lee,  on 
the  other  hand,  gives  empirical  expressions  for  ionization  and 
dissociative  attachment  rates  and  mobilities  that  depend  on 
electron  temperature,  pressure  and  pressure  respectively.  Using 
electron  temperatures  similar  to  Thompson,  and  pressures  of 
about  0.4  torr,  we  obtain  ionization  and  dissociative  attachment 
rates  on  the  order  of  1 0*”  and  10^  .  Thompson  indicates  that 
his  measurements  were  taken  at  0.04  torr,  so  one  would  expect 
smaller  ionization  and  dissociative  attachment  rates,  however, 
since  both  are  directly  proportional  to  the  pressure,  one  would 
not  expect  numbers  as  small  as  are  indicated  by  the  code. 

Likewise,  there  is  a  considerable  difference  in  the 
mobilities  used  in  each  paper.  Lee's  are  about  three  orders  of 
magnitude  larger  than  Thompson's.  Again,  this  seems 
inconsistant  in  that  one  would  expect  smaller  mobilities  with 
Lee  since  he  uses  larger  pressures  than  Thompson. 

If,  however,  the  appropriate  changes  are  made  in  Lee  and 
Thompson,  the  profiles  produced  by  the  code  for  Lee  begin  to 
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take  the  appearance  of  the  profiles  produced  by  the  code  for 
Thompson.  The  changes  necessary  to  do  this  are:  1.  Lee's 
equations  must  be  changed  from  cylindrical  coordinates  to 
Cartesian  coordinates.  2.  Lee's  mobilities  and  ionization  and 
dissociative  attachment  rates  must  be  used  in  Thompson.  3. 
Associative  detachment  must  be  ignored  in  Lee.  4.  The  same 
value  of  h  must  be  used  in  both  codes.  These  changes  put  both 
codes  in  cartesian  coordinates  and  start  both  off  with  the  same 
set  of  parameters.  The  difference  between  the  profiles  found 
from  the  Lee  and  Thompson  codes  decreases  as  h  increases.  At 
h«40,  the  difference  between  both  the  electron  and  negative  ion 
profiles  found  from  the  two  codes  was  less  than  about  3% 
throughout  the  normalized  tube  radius.  Pigs.  7  and  8  show  the 
results  for  h*10  and  40.  Profiles  do  not  extend  to  the  wall 
(z*l)  because  no  variation  of  parameters  was  performed  since  the 
only  purpose  of  the  comparison  was  to  see  whether  the  two  codes 
gave  similar  profiles  given  the  same  initial  parameters. 

Without  associative  detachment/  then/  Lee  and  Thompson  seem  to 
agree  reasonably  well,  however/  keeping  associative  detachment 
in  Lee  and  including  it  in  Thompson  does  not  give  similar 
profiles. 

Boundary  Conditions 

Boundary  conditions  play  an  important  role  in  determining 


the  exact  behavior  of  the  charged  particle  profiles.  Requiring 
the  fluxes/  electric  fields  and  gradients  of  the  charged 
particles  to  be  zero  on  axis  (z*0)  seems  to  be  reasonable  when 


one  considers  the  symmetry  of  the  problem.  Ingold  and  Lee, 
however,  further  constrain  the  problem  by  requiring  the  particle 
densities  to  go  to  zero  at  the  wall  with  the  assumption  of 
quasineutrality  valid  over  the  entire  region.  Another  such 
boundary  condition  is  certainly  necessary  to  adequately  solve 
the  problem,  however,  there  has  been  debate  as  to  whether  the 
condition  of  zero  charge  density  and  quasineutrality  at  the  wall 
is  appropriate. 

Forest  and  Franklin  (Ref.  3)  and  Edgley  and  von  Engle  (Ref. 
2)  argue  that  a  sheath  region  must  exist  at  the  wall  where 
quasineutrality  cannot  be  obeyed.  The  reason  for  this  is  that, 
in  the  sheath  region,  electron  diffusion  is  governed  by  the 
electrons  thermal  velocity  which  is  much  greater  than  both  the 
ion  thermal  velocity  and  the  ion  drift  velocity.  Thus,  once 
inside  the  sheath  region,  electrons  will  diffuse  to  the  wall 
much  faster  than  positive  ions.  So,  to  obtain  equal  fluxes  at 
the  wall,  there  must  be  more  positive  ions  in  the  sheath  region 
than  electrons.  This  arguement  seems  to  preclude  the  existance 
of  negative  ions,  which  is  understandable  for  Forest  and 
Franklin  since  they  only  consider  the  two  component  case.  The 
reason  is  unclear,  though,  for  Edgley  and  von  Engle.  Their 
results  show  that  the  negative  ions  migrate  toward  the  axis  and 
very  few  make  it  to  the  wall,  but  use  of  this  knowledge  in 
setting  up  the  initial  boundary  condition  would  be  an 
inappropriate  approach  to  the  problem. 

For  the  sake  of  arguement,  we  shall  assume  the  approximation 


is  valid  and  continue.  Thus,  at  the  wall,  Edgley  and  von  Engle 
can  equate  the  electron  drift  velocity  to  the  electron's  thermal 
velocity  since  the  electrons  are  close  enough  to  reach  the  wall 
through  the  random  motion  of  thermal  velocity.  Using  simple 
kinetic  theory,  Edgley  and  von  Engle  write  (2:379) 


where  Cj  ^positive  ion  sound  speed,  ‘V'c  *  electron  drift 
velocity,  y«  A- /n~«  ,  Te»  electron  temperature  and  m»  mass. 

Edgley  and  von  Engle  perform  an  integration  routine  similar 
to  Lee  except  they  do  not  assume  quasineutrality  and,  thus,  have 
seven  equations  in  seven  unknowns.  Equation  4-1  provides  a 
boundary  condition  on  the  normalized  electron  velocity  at  the 
wall.  Integration  of  the  equations  proceeds  to  the  v  i.1  where 
the  normalized  electron  velocity  must  be  equal  to  the  ratio  of 
the  particle  masses  as  given  in  4-1.  If  this  condition  is  not 
satisfied,  parameters  are  varied,  as  was  done  for  Lee,  until  the 
boundary  condition  is  met. 

One  possible  criticism  to  this  approach  is  the  use  of  4-1  as 
a  boundary  condition  on  the  equations  governing  the  physics 
throughout  the  rest  of  the  plasma.  Outside  of  the  plasma 
sheath,  the  drift  velocity  of  electrons  is  governed  by  outward 
diffusion  and  an  electric  field  created  by  a  small  deviation 
from  strict  charge  neutrality.  Near  the  wall,  though,  the  drift 
velocity  out  of  the  plasma  is  governed  by  the  electron's  thermal 


velocity.  Thus,  the  equations  which  govern  the  physics  of  the 
sheath  are  different  than  those  which  apply  throughout  the  rest 
of  the  plasma,  yet  von  Engle  and  Edgley  use  the  sheath  condition 
at  the  wall  as  a  boundary  condition  on  the  equations  used  in  the 
integration  routine. 

A  question  also  arises  as  to  whether  the  plasma  sheath  is 
worthy  of  consideration  at  all.  Typically,  the  thickness  of  a 
plasma  sheath  is  on  the  order  of  a  Debye  length,  where,  beyond 
this  distance  the  plasma  remains  undisturbed.  From  Ref.  8  we 
have  the  following  expression  (8:71) 

\x  r  \l GoWr/**  a*  (4-2) 

where  Xjl  is  the  Debye  length  and  He  is  the  electron  number 
density.  Using  Lee's  numbers,  we  find  that  the  number  of  oxygen 

if 

molecules  per  cubic  centimeter  is  on  the  order  of  /0 
Assuming  the  ratio  of  electrons  to  oxygen  molecules  to  be  no 
less  than  1  to  100000,  we  obtain  a  maximum  value  for  Aa  of 
about  0.01  cm.  This  represents  only  1/100  the  radius  of  a  1  cm 
tube. 

Conclusion 

Three  main  conclusions  can  be  drawn  from  this  presentation. 
First,  if  the  same  parameters  and  coordinate  systems  are  used  in 
both  Lee  and  Thompson,  and  associative  detachment  is  ignored, 
then  Lee's  profiles  looked  like  Thompson's;  but,  when 
associative  detachment  is  included  in  both,  the  two  are  not 


similar.  This  may  indicate  that  Thompson's  relation  is  valid 
only  in  limiting  circumstances  (e.g.  when  associative  detachment 
can  be  ignored).  Thompson's  experimental  measurements  are  taken 
at  about  1/10  of  the  pressures  Lee  uses  for  his  profiles  and, 
since  Lee  states  that  the  associative  detachment  rate  is 
proportional  to  the  oxygen  density  and  negative  ion  density,  it 
is  possible  that  associative  detachment  was  not  a  major  factor 
in  Thompson's  experiments.  Some  indication  as  to  where  Thompson 
and  Lee  begin  to  predict  significantly  different  results  can  be 
found  by  comparing  the  two  codes  when  associative  detachment  is 
included  in  Lee  but  not  in  Thompson.  This  was  done  for 
different  values  of  pressure  and  h  and  the  results  are  given  in 
the  table  on  the  following  page.  The  difference  between  the 
normalized  electron  densities  c.nd  the  normalized  negative  ion 
densities  predicted  by  each  code  were  evaluated  at  each  of  the 
100  mesh  points  (except  the  first  and  last) .  The  maximum 
difference  in  the  normalized  electron  densities  is  given  by 
MAXX.  MAXY  represents  the  maximum  difference  in  the  negative 
ion  densities.  DX  and  DY  represent  the  average  difference  in 
the  densities.  X  and  Y  are  the  normalized  electron  and  negative 
ion  density,  respectively,  at  the  wall,  as  given  by  Thompson's 
code.  Values  of  h  larger  than  about  10  resulted  in  rising 
negative  ion  profiles  as  calculculated  in  the  Lee  code. 

Secondly,  the  fact  that  the  power  series  approximation  used 
by  Lee  near  z»0  gave  excellent  agreement  with  the  integration 
routine  out  to  greater  than  z=0.9  indicates  that  numerical 


COMPARATIVE  DIVERGENCE  OP  THOMPSON  WITH  LEE 


H 

P(torr) 

DX{10-4) 

DY(10"4) 

MAXX ( 10  2 ) 

MAXY ( 10-2) 

X 

2 

.02 

.831 

18.2 

.185 

3.86 

.99 

2 

.04 

2.62 

66.7 

.482 

13.54 

.96 

2 

.06 

3.16 

116. 

1.12 

18.37 

.85 

2 

.08 

19.6 

117. 

5.12 

18.8 

.59 

5 

.02 

1.08 

18.1 

.238 

3.92 

.99 

5 

.04 

4.00 

71.3 

.843 

15.4 

.98 

5 

.06 

7.35 

158. 

1.28 

33.0 

.95 

8 

.02 

1.19 

17.6 

.261 

3.84 

1 

8 

.04 

4.52 

70.4 

.970 

15.33 

.99 

8 

.06 

9.20 

157. 

1.85 

34.2 

.97 

.96 


w 


integration  of  equations  C-l  probably  isn't  necessary.  The 
Taylor  series  expansion  used  near  the  wall  could  also  be  used  in 
conjunction  with  the  power  series  since  the  power  series  begins 
to  diverge  near  the  wall.  A  possible  problem  may  occur  in  that 
Lee  only  keeps  the  first  term  of  the  Taylor  series  expansion, 
thus  severely  limiting  the  distance  from  the  wall  where  the 
approximation  is  valid.  This  can  be  overcome  by  writing  the 
expansions,  C-24,  as  infinite  series  and  developing  recursion 
relations  similar  to  C-ll,  C-14,  C-16  and  C-23.  This  was  done, 
but  a  converging  series  could  not  be  obtained.  This  may  be  due 
to  a  simple  algebraic  error,  but  time  limitations  prevented 
further  work  on  this  area. 

Finally,  we  note  the  numerous  problems  that  exist  with 
Ingold's  work.  The  fact  that  B-5  can  not  be  properly 
established  casts  doubt  on  his  entire  development  since  it  is  a 
cornerstone  of  his  work.  Also,  it  was  shown  that  a  completely 
different  solution  could  be  obtained  without  constraining  the 
problem  any  more  than  does  Ingold.  This  seems  to  indicate  that 
Ingold's  solution  may  be  somewhat  arbitrary.  Also,  as  was 
pointed  out  earlier,  the  numerical  solution  developed  for  Ingold 
did  not  give  profiles  similar  to  Ingold's.  This  may  indicate  an 
error  was  made  in  the  development  of  the  numerical  solution,  but 
the  fact  that  the  code  gave  the  proper  results  in  the  two 
component  limit  (as  did  Ingold's  analytic  equations)  serves  as 
an  argument  against  that.  In  fact,  all  numerical  solutions  and 
analytic  equations  reduced  properly  in  the  two  component  limit 


except  Thompson's  relation  for  the  electric  field  (3-30) . 


Recommendations 

Many  areas  of  this  study  could  benefit  from  further  work. 
Among  these  is  the  re-developed  solution  to  Ingold,  done  in 
chapter  2,  using  the  relation 

+by,)-y"(c.x'+Ayl) 

If  another  boundary  condition  at  the  wall  could  be  found, 
then  numerical  profiles  could  be  developed  in  a  manner  similar 
to  the  method  employed  in  appendix  F.  Also  the  alternate 
solution,  2-32,  cold  be  solved  analytically  to  give  charged 
particle  profiles. 

Further  work  could  also  be  done  to  finish  the  developement 
of  an  infinite  Taylor  series  expansion  at  the  wall  as  described 
earlier.  This  could  then  be  coupled  with  Lee's  power  series 
expansions  to  give  a  complete  solution  throughout  the  entire 
discharge  tube,  eliminating  the  need  for  numerical  integration. 
It  may  not  even  be  necessary  to  use  terms  higher  than  first 
order  in  the  Taylor  series  expansions  since  the  power  series 
expansions  were  seen  to  converge  out  to  z=.95,  but  this 
assumption  should  not  be  made  without  further  investigation. 

An  important  final  recommendation  is  to  compare  the  various 
approaches  with  actual  experimental  data  since  a  theory  isn't 
worth  much  if  it  falls  short  of  reality. 


1.  Abramowitz,  Milton  and  Irene  A.  Stegun.  Handbook  of 
Mathematical  Functions.  New  York:  Dover  Publications  Inc., 
1972. 

2.  Edgley,  P.D.  and  A.  von  Engel.  "Theory  of  positive  columns 
in  electronegative  gases,"  Erae^Royal  Society  London.  Hi: 
375-387  (1980). 


3.  Forest,  J.R.  and  R.N.  Franklin.  "The  theory  of  the  positive 

column  including  space-charge  effects,"  British  Journal  of 
Applied  Physics,  1:  1357-1368  (1968) 

4.  Holt,  E.H.  and  R.  E.  Haskell.  Foundations  of  Plasma 
Dynamics.  New  York:  the  Macmillan  Company,  1965. 

5.  Hubin,  Wilbert  N.  Basic  Programming  For  Scientists  and 
Engineers.  Englewood  Cliffs,  New  Jersey:  Prentice-Hall, 
Inc.,  1978. 

6.  Ingold,  J.  "Schottky  Theory  of  the  positive  column  with 
electron  attachment,"  33r  Gaseous  Electronics  Conference, 
(1980) 

7.  Lee,  D.A.  and  Frank  D.  Lewis.  "Multiple  Equilibria  in  the 
oxygen  positive  column,"  Journal  of  Applied  Physics.  51: 
4699-4705  (Sept.  1980) 

8.  Seshardri,  S.R.  Fundamentals  of  Plasma  Physics.  New  York: 
American  Elsevier  Publishing  Co.  Inc.,  1973. 

9.  Thompson,  J.B.  "Electron  energy  distributions  in  plasmas," 
Proc.  Royal  Society.  503-518  (1961) 

10.  Thompson,  J.B.  "Negative  Ions  in  the  Positive  Column  of 

the  Oxygen  Discharge,"  Proc.  Royal  Society.  21:  818-821 

(1959)  . 


APPENDIX  A 


Thompson  begins  with  the  diffusion  equations 


f.  -  “0 

r&  ~  "Oe  £ 


(A-la) 
( A-lb) 
(A-lc) 


which  he  writes  as 


=  -&;  vru 

-  -&*  WU 


•=  -n* 


d;  vne 


(A-2a) 

(A-2b) 

(A-2c) 


by  defining  ambipolar  diffusion  coefficients  as  follows: 


D*=  (N  ro»rvz°<y)  c  i+  >  -1 

1(1  ♦.<*)(  I  «■ 

P>  1±-X  £  - *  1 

L  *  Ac  1  J 

D*-bf  - ^  -V-*  - 1 

^  *L  i  -*■ /a«  J 


(A-3a) 


(A-3b) 


(A-3c) 


To  obtain  A-3,  we  begin  by  stating  the  assumptions  made  by 


Thompson,  namely,  r,=r.  and  |\;»\.+  kic  .  Therefore,  A-la 


be  written  as 


(A-4) 


Subtracting  A-lb  from  A-4  we  have 


le*  f/A.n.E •OJhtTjtA.n+E  (A-5) 

Multplying  A-5  by  and  A-lc  by  ^  and  adding  the 

two,  we  have 

(yu-w.  (A-6) 

Dividing  through  by  the  equation  becomes 

*-^C*-n))?nc  (a-7) 

Rearranging  and  factoring  out  a  -b,  we  have 


Employing  the  Einstein  relation,  1-8,  and  assuming 
the  equation  becomes 


r  *  -^4-^  £  *m)£7 

Ll  +  Si-rCl  ¥~>)/ux4_  t*/*~/uc.  J 


(A-10) 


To  obtain  eq.  A-3c,  we  must  require  the  numerator  in  A-10  to  sum 
to  I  t  }f  l”2«cY  .  This  can  be  done  one  of  two  ways:  the  first 
is  to  assume  yU+  ~M-  which  is  not  always  true.  If  the 
assumption  is  not  made,  we  have 


-  J±Z  r 

A.*.  7 Me 


,  .1 1 

V*le 


A- 


or 


(A-ll) 

With  this  proportinality  relation,  A-10  reduces  to  A-2c. 

To  obtain  A-2b,  we  begin  by  subtracting  A-lc  from  A-4  to  get 


fj=-D+Wfy,fw+E E 


(A-12) 


Multiplying  A-12  by  rt-  and  A-lb  by  Ac»  an<* 

adding  we  have,  after  dividing  by  • 


A* 


-  (  % *- 1  i  *•  £  •<  £>. 


(A-13) 


A* 


/  *  /**+  C  t  +-  -cj  f-  /%- 


where  the  Einstein  relation  and  the  fact  that  f7*U  has 

again  been  used.  Again  we  want  the  numerator  to  reduce  to 


l  .  Thus# 


Y-t-  ^  y 


Vn- 


3  ^  (*r  *  ir'vNi ■*  -  r+i«.r 

7n-  '  /^+ 


=  /  1-  Y  i-  7<*ir 


^Le  W(|f*  £*>/  +  <£!*' 
7n_ 


(A-14) 


At  this  point#  we  can  see  that  assuming  ju-  will  not 
provide  a  solution.  Instead#  we  obtain  the  same  relation 
previously  found#  i.e.^^si  Y 

To  establish  A-2a#  we  add  Al-b  to  A-lc  and  multiply  through 
by  yU+n*  to  get 


^M4.n+  H  *  Dyi*hn+  7h~  E 


(A-15) 
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•  r* 


'•lr 


l  f  < 


Again#  using  the  Einstein  relation  and  the  fact  that  ft 
and  Vn«. -V/L+Vrt*  #  we  obtain 

+  I  K 

Using  the  relation  7n*  6  °s  #  we  get 

P_  K_  fnitockc* £i(<+">nk:i'l**%- 
£■*'  *  n{  /  + /"*u++)*  % it 

Factoring  out  a  (l  f  ^Z)/(l b  **)  •  we  have 

P  -  -h  «7/j  f — 2E»  LLll  ^  y  —  I 

If'  ^  +1  ( t  r<x  ( ^O+^t-w^yA* 


APPENDIX  B 


Derivation  of  Inaold's  equations 


Ingold  begins  with  the  diffusion  equations: 


C- 

-J±-*  E 


Together  with  the  continuity  equations: 

d\ 

-  o  «-*. 

At 

•If-  X  y* 


and  an  assumption  of  quasineutrality , 


Ingold  obtains 


f6,w©s)  +  VJn«-- 

Through  further  manipulation  of  equations  B-] 
Ingold  is  able  to  obtain  the  following  relation: 


(1 

(1 

(1 


(1 

0 
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c 

and  B-2 


(B-5) 


nJ 


v.  /**  /**■-  /**.  / 


Ingold  fails  to  give  any  details  as  to  how  this  last 
equation  was  obtained  however.  Dr.  Alan  Garscadden  and  Lt.  Col, 
William  Bailey  offered  some  insight  by  developing  a  similar 
relation.  The  derivation  is  as  follows. 

Rearranging  B-lb  we  have 


n 

JvG<ne+*-)E  =  -*•*&*/ 

7 


(B-6) 

Rearranging  and  adding  B-la  and  B-lb  and  replacing  the  electric 
field  term  with  B-6  we  have 


^  +  £:  W')i-  &  He'  =  O 

/*v  a.«  ^  3 


5 


(B-7) 


Differentiating  this  last  equation  and  substituting  into  B-2 
gives  us 


(S  »s h*" 


(B-8) 


Assuming  we  get 


(&e^5)  -£H*  <0 


(B-9) 


which  verifies  B-4 


Differentiating  B-2b  and  making  the  proper  substitutions  we 


have 


(B-10) 


Differentiating  B-la  and  B-lc  and  adding  both  to  B-10  we  obtain 


(B— 11) 


Rear ranging ,  we  have 

/si  A  - 


Ea 


T.  M- /  (  I  f  j 


(B-12) 


Differentiating  B-lc  and  substituting  in  B-12  we  obtain 


0  K  " (iKl+K* 


(B-13) 


This  differs  from  Ingold  in  that  it  relates  the  ratio  of  second 
derivatives  to  first  derivatives  whereas  Ingold  relates  the 
ratio  of  first  derivatives  to  undifferentiated  terms.  It  can  be 
shown  that  B-16  reduces  to  B-5  if  the  ratio  can  be  written  as 
two  separate  equations,  i.e. 


0* 

t  /*<■  ly*"*  M~ 


but,  as  was  shown  in  the  derivation  of  B-16,  this  can  not  be 
done  since  the  term  common  to  both,  is 

cancelled  out  when  the  ratio  is  taken.  Thus,  B-5  cannot  be 
properly  established  and,  since  it  is  an  essential  relation  to 
Ingold's  development  (as  will  be  shown  later),  his  final  results 
must  be  viewed  with  skepticism. 

Let  us  now  begin  the  task  of  tracing  out  Ingold's 
development.  It  should  be  noted  that  only  those  equations 
specifically  referred  to  as  Ingold's  were  obtained  from  Ingold's 
notes.  All  other  equations  are  conjecture  as  to  how  Ingold's 
equations  were  arrived  at. 


He  begin  by  integrating  B-2  with  respect  to  z 


IX-J-X)  N«. 

£_  --  \  Me 


Thus,  B-la  can  be  written  as 
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x  -o  /.tAVc.x 
Similarly,  B-lc  can  be  written  as 


^  It- 
C  1  ■»  ^ 


Equating  B-18  to  B-19, 


&e*i  * 


we  obtain  Ingold's  relation 

v ~x_  xv  . /^  xVa  xy _ 

-0  /‘•AfX'jry 


Substituting  B-17  into  B-7,  we  have 


£*+ ^  <.e,<4Ve,n.Veene  =0 

**+  ■>»  #«,  /A*  /*“ 


/*+  *0  #«,  /A*  /* 
Solving  for  XV  ,  we  obtain 


p  -C©4^)/i«*x  i-ie5n-,y' 
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(B-17a) 

(B-17b) 

(b-17c) 


(B-18) 


(B-19) 


(B-20) 


(B-21) 


(B-22) 


Substituting  B-22  into  B-20,  we  get  Ingold's  next  relation 


Differentiating  B-23  and  evaluating  it  at  z*0  (where  A>ay^'sl 
and  ^sry##*0  ),  we  obtain  the  following  Ingold  equation 
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:'-± 
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Defining  a,b,c  and  d  as 
©e  /  d«  \  X 


B-5  can  be  written  as 


*'Un- 


(B-; 


(b-: 


Differentiating  B-25  and  evaluating  it  at  z»0  we  have 


&''(>*•  k)  -  (B-: 


It  will  be  shown  later  that 
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67 


Thus,  we  obtain  Ingold's  relation 
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(B-29) 


But,  as  was  pointed  out  earlier,  B-5  can  not  be  properly 
established.  Instead,  B-16  must  be  used  which  yields 


(B-30) 


instead  of  B-5.  Since  B-29  is  needed  in  the  next  step  of 
Ingold's  development  and  B-26  in  derivations  following  that,  all 
results  past  this  point  must  be  considered  dubious  at  best. 

Assuming  B-29  is  true,  and  substituting  it  into  B-24  and 
simplifying,  we  obtain  one  of  Ingold's  key  results 


•O-X  X  _ 

/*«■ 


(B-31) 


Let  us  now  establish  the  following  two  equations  given  by 


C+d-  I  +■ 


M-t.  +  K -o 

i  • 


(B-33) 


From  B-25  and  B-31,  we  have 


^  + &*-] 


(B-34) 


We  now  factor  out  an  A  which  will  cancel  an  A  that  is 
factored  out  of  the  expression  c+d,  i.e.  we  essentially  redefine 
a,b,c  and  d  such  that 


\  Gon-b  y3 


(B-35) 


Therefore,  B-34  reduces  to 


(B-36) 


Placing  all  terms  over  a  common  denominator  of  yu*yU._ 

factoring  out  a  A«- +■/*-* -0  f  we  have 

/++/+- 
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Again  factoring  out we  obtain 


(B-39) 


Again  we  redefine  a,b,c  and  d  such  that  the  terms 
cancel  in  equations  B-37  and  B-38. 

We  will  now  show  that 


U/jt  Gi?\  M*/±r  /«+ *««  +/+*'ie*y«y.y\-'> tA - *!-» 

'  )Vo  /  Ac  i  K>« 

/»«  H*,  ^-o 


(B-40) 


thus  proving  that  B-37  and  B-39  reduce  to  B-32  and  B-33. 


Since 


^  At +/ju.- 

££.  4.  A«.he»  yu-A-o 

/.  -  ttcua  /4.  Ac»  A-^e» 


*  _  t^l  <Lj 

/*<*<* '  M^TyU-n-C 

yU~  A*.0 


we  can  write 

f  1  .  ^-qN  ^Z ±£L  - 

V  ^  Aetfj  +  *S=*  ~  .  ft  4. ,  *  . 

'  >u-  Atn«  fA- n-* 


(B-41) 


(B-42) 


Therefore 


a»/a- 

4-  h*® 
A- 


(B-43) 


which  verifies  B-40.  Although  keeping  the  terms  we  cancelled 
out  by  redefining  a,b,c  and  d  would  have  given  expressions 
different  than  B-31  and  B-32,  the  important  relation  needed  fo 
Ingold's  development,  equation  B-28,  would  still  be  valid. 

Continuing  Ingold's  development.  Ingold  next  assumes  that 


LA.  +- 


(B-44) 


where  n  is  some  function  of  the  spatial  coordinate  and  f  is  a 
constant.  It  follows  then  that 

^ ia  t  -f  -  ky  _  (^U)x  4-bg  (B_45 

C kyc  C-jc  1-oIm. 


A  u  _  ~  ^  (tvcl4)J  yi  +  (lo~ X$)  U 

(k*  +•  ci.  VA 


Ingold  then  makes  the  seemingly  arbitrary  assumption 


--f(  tt-el-f)  -0 


(B-47) 


which  £ixes  the  value  of  f  at 


(B-48) 


This  allows  Ingold  to  write  B-46  as 
dii  _  V>s  M  _ 

(jlx  Ax  *  f  k 


(B-49) 


where  m=b-df,  n*c+df  and  p®d.  The  solution  to  B-49  is 


+■ 


je- 

»\-n 


U 


(B-50) 


Since  y»x*l  at  z*0,  we  can  see  from  B-44  that  u*l-f  at  z*0. 
Using  these  values  for  x,y  and  u  determines  A,  which  yields 
Ingold's 


i-7trn0-() 

(i  -f)n/m 


(B-51) 


Substituting  B-50  into  B-44  we  obtain  Ingold's  expression  for  y 


'y  -  Cu  +•  D  4 

where 

CM-f  •  £> -Bf+I  ;  & 

From  B-9,  using  the  variables  x  and  y,  we  have 

where 

,  i_  +  — 

K  ~  A*  /*»  /*- 

Thus,  taking  the  second  derivatives  of  B-50  and  B-52  and 
substituting  into  B-54,  we  obtain  Ingold's  key  result: 

((A+l^su'V  d]u\(a*&C>«-»)u 


APPENDIX  C 


Derivation  of  Lee's  Recursion  Relations 


Lee's  four  first  order  differential  equations  including 
boundary  conditions  are  (7:4701) 


(C-la) 

i'Liyt)'-  ),•<$/>.)  y> 

(Olb) 

y.'+^y/ 

(Ole) 

>'y.  -  cy.yj 

(C-ld) 

and 

y,«))  =  0;'yilp)=0^  y^O-0 

(02) 

where 

v 

1*  “  /l«o  /  7l  “  fl«o 

V  v  - 

(03) 

and 

M  t)+  \  JO 

T  *  V©e  ;  *  =  - 

****  e.  ; c-  -^r 

(04) 

Also,  the  derivatives  are  in  terms  of  the  normalized  coordinate, 
z,  and  j  represents  the  charged  species  flux  density.  The 
momentum  transfer  frequencies  are  given  by  Lee  as  (7:4703) 

sec'  ** r '•  xP 


(C-5) 


Since  Lee  uses  cylindrical  coordinates,  a  singularity  is 
encountered  at  z»0.  To  get  around  this,  Lee  uses  power  series 


approximations  near  zero  given  by  (7:4703) 


where 


V. :  ;  y.=  ‘j.tS-t.,*1'1' 

9*  ~  Xk-I 

)3~ 


(C-6) 


(07) 


Substituting  03  into  Ola,  we  have 


£(?  2lk’V 


Redefining  the  dummy  variable  of  summation  on  the  right  side  as 
k*k-l,  we  obtain 


K-, 


*  +z>; 


(09) 


If  we  include  po  and  C j,  in  the  summation  on  the  right  side 
of  09,  we  have 


Again  redefining  t-^e  dummy  variable  k*  Jt-*'/  and  equating  like 
powers  of  z,  we  obtain  Lee's  eq.  4.5  (7:4703) 

i* e.i.i-  <c-u> 

Substituting  C-3  into  C-lb,  we  have 

i(?su.yM+  (c-i2) 

where  p0  and  <£0  have  again  been  included  in  the  summation 
of  y t  and  .  We  redefine  k  =  k-*-l  on  the  left  side  of 

0-12  to  obtain 

go  ..  «*•  *1.  *5  Zk  (PP  T-k 

1 W  *  ?**.,(.«■♦  (c-13> 

Equating  like  powers  of  z,  we  obtain  Lee's  eq.  4.6  (7:4703) 

~  W)  ‘ '  *  *  *  *  *  (C-14) 

Substituting  C-3  into  C-lc,  we  have 

It  t2(r  4U_I 


Equating  like  powers  of  z#  we  obtain 


P«<  ~  ~  ii,  **”  ^  ) 


(C-16) 


which  is  identical  to  Lee's  eq.  4.7  (7:4703). 
Finally#  substituting  C-3  into  C-ld#  we  have 


(C-17) 


Multiplying  the  first  two  terms  together#  we  obtain 


(C-18) 


where  we  have  changed  the  summation  variable  from  k  to  d  in 
the  second  term  for  clarification.  Simplifying  the  second  term 
in  C-16,  we  have 


ck$  n  -? 00 


(C-19) 


where  .  Since  we  are  interested  in  equating  like  terms 

of  z,  we  want  all  combinations  of  that  yield  the  same  m 

All  such  combinations  can  be  expressed  as 


for  any  given  m.  Likewise, 


becomes  -C  2l  C2-kJ  ^ 

becomes  .-V*  (C-21) 


Therefore,  upon  substituting  m  for  k  in  the  remaining  sums  and 
equating  like  coefficients  of  z,  we  obtain 


Since  and  yo«x|  we  have,  after  changing  summation 

variables  from  m  to  n  and  from  k  to  1  and  simplifying, 

~  t11*-*-) )  ]  (c'23) 

This  differs  from  Lee’s  eq.  4.8  (7:4703)  in  that  the  third  term 
on  the  right  side  of  the  equation  contains  an  extra  £  ,  i.e. 

/i-t  .  . 


After  discussing  the  discrepancy  with  Or.  Lee,  it  was  discovered 
that  the  missing  JL  was  merely  a  typographical  error  in  his 
paper . 

The  singularity  at  the  wall  is  dealt  with  by  the  formal 


approximation  (7:4703) 


y.-a.ft-'a)  f  oCl-^)x 

(C-24a) 

(C-24b) 

yj=c.  +c((i-ar)  +  o(i-*jx 

(C-24c) 

yH  -d*  +  cA.O-3) 

(C-24d) 

Lee  then  substitutes  C-24  into  C-l  to  obtain  (7:4703) 

\jjic.fx-t)  yH:cL  Cz-t) 

with 

Co  -  D^Cb-aXmM^^.&o-ct)*] 

^  -  5“[a  slo(  (l  -O')  -  ct,6i, + 2-<rl>,)3 

D  =  A.-ACe-A)  -t-BCl, 

where  terms  of  order  (I  -  are  neglected.  These  equations  will 
now  be  derived  for  verification. 

Substituting  C-24  into  C-lc  gives 


(C-25) 


(C-26a) 

(C-26b) 

(C-26c) 


(c-27) 


simplifying,  we  have 


(C-28) 


Now  substituting  024  into  Old,  we  obtain 

-  + r  V  A  ( I  -  -  (13 (l-i 

Substituting  028  into  029,  we  have 

-A,k,Cl  -*)+ 

Neglecting  terms  of  order  and  solving  for 

obtain 

_  A,  C8-AJ(f<,  -t-lrla,)  i-<tAi  BO  -<D 
«.A(8'A)  +BCk, 


Thus  verifying  026a.  Now  solving  029  for  C*+  C,0 
substituting  it  into  028,  we  get 


JLf«l  (i-a)-  A[(BA)yCi-t)[«U'*«l.Ci-x>V‘M.(i 

1  B 


-1 

(029) 

*»]/s 

(030) 

Cm  ,  we 

(031) 

b )  and 

jj] 

(032) 


Multiplying  through  by  ,  neglecting  terms  of 

order  and  solving  for  d*  ,  we  obtain 


a,A£0-a)  +■  BCl», 


(C-33) 


Thus  verifying  C-26b. 

We  must  now  attempt  to  verify  C-25.  We  see  that  the 
expressions  for  \jt  and  yx  in  C-25  are  the  same  as  C-24  but 
those  for  yi  and  ^  change.  Let's  begin  with  C-lb. 


ZYl'*'/'.  ■**(>'  f  >0 

Adding  and  )x  to  both  sides,  we  can  write 

-  Y,(i-y  i-'/,  - 'f,  Y  v=  f  Y')  6-W'Ym' 


(C-34) 


(C-35) 


Substituting  C-24  into  C-35  and  ignoring  terms  of  order  , 


we  obtain 


(C— 36) 


(iXr  x  06-*)^<W< 


( C— 37 ) 


Since  the  left  side  depends  on  (l~?<)  and  the  eight  side  does 
not,  we  must  require 


(C-38) 


It  is  important  to  note,  at  this  point,  that  the  sum  of  the 
coefficients  on  the  left  side  cannot  be  equated  to  zero  since, 
strictly  speaking,  equations  024  should  be  expressed  as  an 
infinite  series  in  powers  of  (l-a)  .  Thus,  the  derivation 

of  y4  would  be 


I  "f'cp/ivj 


(039) 


where  el*,  would  have  to  be  included  on  the  left  side  of  037. 
(Credit  for  the  above  information  is  attributed  to  Dr.  Lee  for 
his  informative  discussion  on  the  apparent  contradictions 
created  by  ignoring  second  order  terms.) 

Looking  now  at  Ola,  we  have 

^  'ii  +7)  “  ^-X)A]  ,+  ^  *  (040) 

Adding  -  y3'  and  "^((V-Xy)Jy,  to  both  sides,  we  obtain 


(041) 


Substituting  024  into  041  and  ignoring  terms  of  order  ( l—  £) 
we  obtain,  after  rearranging 


[c,  -X)AK  ~  f  -LrCo  (042 

Again,  equating  zeroth  order  terms  in  ,  we  have 

C,  -  Ce  (043 

Using  038  and  043,  we  see  that  O24o  and  024d  reduce  to 

Cl'  2?)  (044 


thus  verifying  025. 


APPENDIX  D 


Explanation  of  Lee's  Code 


List  of  Variables 

A,B,C  -constants  used  in  Lee’s  empirical  formula  for 
calculating  ionization  rate 

Al-4  -constants  used  in  calculating  polynomial 
approximation  for  erfc(Cfc/j^) 

AA,BB,CC  -A/B  and  C  as  given  in  eqs.  C-4 

AAl , BB1 ,  -  b  ,  c  #  d  and  D  as  referred  to  in  024,025  and 

CCl , DD1 ,  C-26 

DD 

ALPHA  -dissociative  attachment  rate  (sec-1) 

BETA  -ionization  rate  (sec-1) 

Cl-11  -constants  used  in  calculating  Lee's  empirical  formula 
for  dissociative  attachment  rate 

DIPl  -difference  between  Y3#  calculated  using  boundary 

condition  equations  and  the  value  calculated  at  the 
boundary  using  the  integration  routine 

DIP2  -difference  between  Y4#  calculated  using  boundary 

condition  equations  and  the  value  calculated  at  the 
boundary  using  the  integration  routine 

DT*  -mesh  interval,  6%  ,  used  in  integration  routine 

EE  -electric  field  calculated  from  electron  diffusion 

equation 

EN  -electric  field  calculated  from  negative  ion  diffusion 

equation 

EP2  -distance  from  wall  where  integration  routine  is  forced 

to  meet  boundary  condition. 

El  -constant  used  in  Lee's  empirical  formula  for 

ionization  rate 

FNERF  -Polynomial  approximation  for  erfc(C,/7i)  (1:299) 
-associative  detachment  rate  (sec-1) 


GAMMA 


H 

Kl#-4# 

I 

K 

L% 

MNUE , 
MNUP,MNUN 

NK 

NN 

N02 

P 

P#(L%) 

Q#U%) 

R#(L%) 

S#(L%) 

SIGMA 

SUMl#-4# 

TO 

Tl 

TE 

TN 

Yl#-4# 

Y1D4-4D# 

Y10I-40# 


-on  axis  ratio  of  negative  ions  to  electrons 

-first  derivative  approximations  used  in  integration 
routine 

-do  loop  counter 
-do  loop  counter 

-number  of  terms  retained  in  polynomial  approximations 
given  by  C-6 

-constants  given  by  C-5 
-do  loop  counter 

-number  of  mesh  points  used  in  integration  routine 

-3 

-number  of  oxygen  molecules  per  cm 
-pressure  in  torr 

-coefficients  for  polynomial  approximation  of  Yl# 

-coefficients  for  polynomial  approximation  of  Y2# 

-coefficients  for  polynomial  approximation  of  Y3# 

-coefficients  for  polynomial  approximation  of  Y4# 

-ratio  of  backround  gas  temperature  to  electron 
temperature 

-summations  indicated  in  C-22  and  C-23 
-backround  gas  temperature  (ev) 

-initial  value  of  z  where  integration  begins 
-electron  temperature  (ev) 

-value  of  z  where  integration  ends.  EP2=1-TN 
-  y,  ,  and  as  indicated  by  C-3 

-derivatives  of  Yll-4#  in  power  series  approximations 
-values  of  Yl#-4#  from  previous  mesh  point 


*.  *- \  J V 
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2/ 

NL 


Y13I-43#  -intermediate  approximations  for  Yl#-4#  contained 
within  the  integration  routine 

YPR1#  -evaluation  of  Y1PRIME#  for  calculation  of  Y2# 

YPRIMEl#-4#-derivatives  of  Yl#-4f  in  integration  routine 

-value  of  z  at  which  power  series  approximations  are 
made .  Z«T1 


Z 


Evaluation  of  coefficients  for  power 


series  expansions  of  y^r  y2r  y3»  y4 
using  recursion  relations  Oil,  C-14 ,  C-16 ,  023 


Evaluation  of  y^,  y2 ,  y3  and  y4  using 
power  series  expansions 


Beginning  of  fourth  order  Runge-Kutta  integration 
routine  using  values  of  y^,  y2 •  y3  and  y4 
from  power  series  approximations  as  initial  values. 
Evaluation  at  mesh  point  NK. 


1000  DIM  R# (50) ,S#(50) ,P#(50) ,Q#(50) 

1010  P-.4 

1020  MNUN«4.8693E-16*P 
1030  MNUP«8.8508E-16*P 
1040  MNUE«2.9127E-18*P 
1050  TO*. 15 
1060  TE-2.1 
1070  SIGMA-TO/TE 

1080  N02* (P*1333 ! ) / (8 .313 6*T0* 11 587*1 .66 E-24*100 *100000 1 ) 

1090  C1-1.84E-11:  02*1.15:  03*47.7:  04*181:  05-251:  C6-159: 
07*110 

1100  08*52.6:  09*91.1:  010*64.2:  Cll-47.9:  A-2.68E-09: 
B*4.61E-09 

1110  C*.045:  Al*. 278393:  A2*. 230389:  A3*. 000972:  A4=. 078108 

1120  DEF  FNERF(X)  *1-1/ (l+Al*X+A2*X*2+A3*X*3+A4*X*4)  *4 

1130  ALPHA*N02*Cl*EXP{-3 .8/TEJ/TE* (3/2) * ( .  8862*EXP( (C2/TE) *2) 

* ( 1-FNERF ( C2/TE ) ) * ( 03  -  C4/TE  +  C5/TE“2  -  C6/TE*3 
+  C7/TE*4)  +  08  -  C9/TE  +  C10/TE*2  -  Cll/TE~3) 

1140  EI-12.063 

1150  BETA*N02*A*EXP(-EI* (1/TE  -  C) ) /  (TE~ (3/2) ) 

*1/ (1/TE  -  C) *(1/ (1/TE  -  C)  +  El) -B*EXP(-EI/TE) 

* (TE* ( .5)  +  EI*1/TE* .5) 

1160  PRINT  "alpha*"; ALPHA, "beta*"; BETA 
1170  AA«MNUP*ALPHA/(TE*1.6E-12) 

1180  BB-AA  +  MNUN*ALPHA/(TE*1.6E-12) 

1190  CC»MNOE*ALPHA/(TE*1.6E-12) 

1200  GAMMA* .03*NO2*l . IE-09 

1210  PRINT  "aa»";AA,"bb*";BB,"cc»";CC 

1220  HI*100 

1230  REM  DEVALUATION  OF  COEFFICIENTS  FOR  POWER  SERIES 
EXPANSIONS  OF  Yl,  Y2 ,  Y3  AND  Y4** 

1240  H*20 
1250  L%=-1 
1260  P#(0)*1 
1270  Q#(0)*H 

1280  SUM1*«0:  SUM2#*0:  SUM3#*0:  SUM44-0 
1290  FOR  1*0  TO  15 

1300  R# (2*1+1) *( ( (BETA-ALPHA) /ALPHA) *P# (2*1)  +  (GAMMA/ALPHA) 

*Q#(2*I))/(2*(I+1)) 

1310  S# (2*1+1) * (P# (2*1)  -  (GAMMA/ALPHA) *Q# (2*1) )/(2* (1+1) ) 
1320  FOR  J*1  TO  I 

1330  SUM1#»SUM1#  +  S# (2*J-1) *P# (2* (I+l-J) ) 

1340  SUM2#*SUM2#  +  2*J*Q# (2*J) *Pt (2* ( I+l-J) ) 

1350  SUM3#*SUM3#  +  2*J*P# (2*J) *Q# (2* ( I+l-J) ) 

1360  SUM4#*SUM4#  +  R# (2*J-1) *Q# (2* ( I+l-J) ) 

1370  NEXT  J 

1380  P# (2*1+2) *({ (BB-AA) *(S# (2*1+1) +SUM1#)  +  SIGMA*SUM2# 

-  SUM3*  -  CC* (H*R# (2*1+1) +SUM4#) )/ (2* ( 1+1) ) 

-  (AA*R# (2*1+1)  +  BB*S# (2*1+1) )/(4* (1+1) ) ) 

/ ( .5+H) 


1390  Q# (2*1+2) «(-(AA*R# (2*1+1)  +  BB*S#(2*I+1) )/(2*(I+l) ) 

-  P# (2*1+2) )/(2*SIGMA) 

1400  L%»L%+1 

1410  NEXT  I 
1420  71*-P«(0) 

1430  Y2#»Q#(0) 

1440  Y3#-0 
1450  Y44-0 
1460  Z-.01 

1470  REM  ‘‘EVALUATION  OP  INITIAL  VALUES  OF  Yl#  Y2,  Y3  AND  Y4 
USING  POWER  SERIES  EXPANSIONS** 

1480  FOR  K«1  TO  L«+l 

1490  Y1#»Y1#+P# (2*K) *ZA (2*K) 

1500  Y1D#»Y1D#+P# (2*K) *2*K*ZA (2*K-1) 

1510  Y2#-Y2#+Q#(2*K)*ZA(2*K) 

1520  Y2D#«Y2D#+Q# (2*K) *2*X*ZA (2*K-1) 

1530  Y3#»Y3#+R# (2*K-1) *ZA (2*K-1) 

3  M0  Y3D#«Y3D#+R#  (2*K-1)  *  (2*K-1)  *ZA  (2*K-2) 

1550  Y4#-Y4#+S#(2*X-1)*ZA(2*K-1) 

1560  Y4D#-Y4D#+S# (2*K-1) * (2*K-1) *ZA (2*K-2) 

1570  NEXT  K 

1580  REM  * *RUNGA-KUTTA  FOURTH  ORDER  METHOD** 

1590  T1-.01 

1600  TN* .99 

1610  NN-100 

1620  DT#«(TN-T1)/NN 

1630  FOR  NK«1  TO  NN 

1640  Y1T#*Y1# :  Y2T#-Y2#:  Y3T#»Y3#s  Y4T#«Y4# 

1650  T4#-T1+NK*DT# 

1660  T3#«T4f-DT#/2 

1670  T2#-T4#-DT# 

1680  REM  * Evaluation  of  Yl* 

1690  Y10#«Y1# 


1700 


1710 

1720 


1730 

1740 


1750 

1760 


1770 

1780 

1790 

1800 

1810 


Kl#«( (BB-AA) *Yl#*Y4#-CC*Y2#*Y3#-Y10#*(AA*Y3*+BB*Y4#)/2) 

/ (Y2t+Yl#/2) 

Y13#«Yl#+Kl#*DT#/2 

K2  #« ( ( BB-AA) *Y1 3  f *Y4  #-CC* Y2  #  * Y3  *- Yl 3  *  * ( AA* Y3  #+BB*Y4  * ) /2 ) 
/ (Y2#+Yl3#/2) 

Y13#»Yl#+K2**DT#/2 

K3#»( (BB-AA) *Y13#*Y4#-CC*Y2#*Y3#-Y13#*(AA*Y3*+BB*Y4#)/2) 
/ (Y2#+Yl3#/2) 

Y13*-Y1#+K3I*DT# 

K4 #» ( (BB-AA) *Y1 3  **Y4  #-CC* Y2  #  * Y3  4-Y13  * * ( AA* Y3  #+BB*Y4# ) /2 ) 
/(Y2#+Y13#/2) 

Yl #-Yl #+ (K1 #+2 *K2#+2  *K3  #+K4  # ) *DT#/6 
YPRIMEl#-(Kli+K2#*2+2*K3t+K4#)/6 
REM  ‘Evaluation  of  Y2* 

Y20I-Y2* 

YPR1#«( (BB-AA) *Y10#*Y4*-CC*Y2#*Y3#-Y10# 
*(AA*Y3#+BB*Y4#)/2)/(Y2#+Y10#/2) 


90 


1820  K1 #- (-AA‘Y3 #-BB*Y4#-YPRl # ) / ( 2  * . 15/2 . 1 ) 

1830  Y23#»Y2#+Kl#*DT#/2 

1840  YPRl*-( (BB-AA) *Y10#*Y4#-CC*Y23**Y3#-Y10# 

* (AA*Y3#+BB*Y4#)/2) / (Y23#+Yl0#/2) 

1850  K2#» <-AA*Y3  #-BB*Y4#-YPRl # ) / { 2* . 15/2 . 1 ) 

1860  Y23#«Y2#+K2#*DT#/2 

1870  YPR1#«( (BB-AA) *Y10#*Y4#-CC*Y23**Y3#-Y10# 

*(AA*Y3#+BB*Y4#)/2)/(Y23#+Y10*/2) 

1880  K3  #» ( -AA* Y3  #-BB* Y4  #- YPR1 # )  /  ( 2  * . 15/2.1) 

1890  Y23#-Y2#+K3#*DT# 

1900  YPR1#«( (BB-AA) *Y10#*Y4#-CC*Y23#*Y3#-Y10# 

* (AA*Y3#+BB‘Y4#)/2) / (Y23#+Y10#/2) 

1910  K4  #- ( -AA*Y3 #-BB* Y4  #- YPR1 ♦ ) / ( 2  * . 15/2.1) 

1920  Y2#-Y2§+(Kl#+2*K2#+2*K3#+K4#)*DT#/6 

1930  YPRIME2l-(Kl#+K2#*2+2*K3#+K4#)/6 

1940  IF  YPRIHE2#>0  THEN  GOTO  2220 

1950  REM  ‘Evaluation  of  Y3* 

1960  Y30#«Y3# 

1970  Kl#» ( (BETA-ALPHA) /ALPHA) *Y10#+Y20#*GAMMA/ALPHA-Y3#/T2# 
1980  Y33#-Y3#+Kl#*DTf/2 

1990  K2#«( (BETA-ALPHA) /ALPHA) *Y10#+Y20#*GAMMA/ALPHA-Y33#/T3# 

2000  Y33#«Y3#+K2#*DT#/2 

2010  K3#»( (BETA-ALPHA) /ALPHA) *Y10*+Y20#*GAMMA/ALPHA-Y33#/T3# 

2020  Y33#-Y3#+K3#*DT# 

2030  K4#» ( (BETA-ALPHA) /ALPHA) *Yl0t+Y20#*GAMMA/ALPHA-Y33#/T4# 

2040  Y3#«Y3#+(Kl#+2*K2#+2*K3f+K4*) *DT#/6 

2050  YPRIME3#«(Kl#+K2#*2+2*K3#+K4#)/6 

2060  REM  ‘Evaluation  of  Y4* 

2070  K1#»Y10#-Y20t ‘GAMMA/ ALPHA  -  Y4I/T2# 

2080  Y43*«Y4#+Kl#*DT#/2 

2090  R2#*Y10t-Y20# ‘GAMMA/ ALPHA  -  Y43I/T3# 

2100  Y43#*Y4#+K2#‘DT#/2 

2110  K3 #»Y10 #-Y2 0 # ‘GAMMA/ ALPHA  -  Y434/T3# 

2120  Y43#»Y4#+K3#*DT# 

2130  K4#-Y10*-Y20# ‘GAMMA/ ALPHA  -  Y434/T4# 

2140  Y4#«Y4#+(K1#+2*K2#+2*K3#+K4#) *DT#/6 

2150  YPRIME4*- (K1#+K2#*2+2*K3#+K4*) /6 
2160  EE«-ALPHA*Y3#*MNUE/(1.6E-19*Yl#*lE+07) 

-  1.6E-19*TE*YPRIME1*/(Y1#*1.6E-19) 

2170  EN»- ALPHA* Y4  #  ‘MNtJN/  ( 1 . 6  E- 1 9  *  Y2  #  *  1 E+0  7 ) 

-  1.6E-19*T0*YPRIME2#/(Y2#*1.6E-19) 

2180  PRINT  T4#,EE,EN 

2190  PRINT  Yl#,Y2#/H,Y3#,Y4# 

2200  PRINT 
2210  NEXT  NK 
2220  H-H-.001 

2230  LPRINT  "H-";H+.001,,,NK-";NK 
2240  EP2-.01 
2250  AA1»Y1#/EP2 
2260  BB1»Y2#/EP2 


2270  DD«AA1*AA* (BB-AA) +BB*CC*BB1 

2280  CC1- (AAl* (BB-AA) * (AA1+2*SIGMA*BB1) +AA1*BB1*BB* (1-SIGMA) )/DD 
2290  DD1— (AA*AA1*BB1* (1-SIGMA) -CC*BB1* (AA1+2*SIGMA*BB1) ) /DD 
2300  DIF1«Y3#-CC1* (1+EP2) 

2310  DIF2«Y4#-DD1*(1+EP2) 

2320  LPRINT  "difl-" jDIFl, "dif2»" ;DIF2 
2330  GOTO  1250 


ALPHA 


-dissociative  attachment  rate  (sec  ) 

C2  -constant  given  by  eq.  B-54 

DIP  -difference  between  the  left  and  right  sides  of  eq.  3-6 

DT  -mesh  interval 

E  -electric  field  (V/cm) 

GAMMA  0./©5 

HI ,L0  -boundaries  on  x  used  in  the  iterative  subroutine  for 
finding  x 

Kl-4  -estimate  of  first  derivative  of  w  in  integration 
routine 

Ll-4  -estimate  of  second  derivative  of  w  in  integration 

routine 

M  -coefficient  of  X*  in  definition  of  w  (M*  ) 

MUE/MUI /  ^ 

MON  -mobilities  ^ua  yt**  and^u,. 

N  -number  of  mesh  intervals 

NC  -coefficient  of  x  in  definition  of  w  (NOl) 

NO  -ionization  rate  (see"*) 

Q  -final  estimate  of  the  first  detivative  of  w  at  the  end 

of  each  iteration  in  the  inegration  routine 

Q1  -initial  value  of  Q 

R  -ratio  of 

SE,SI,SN  -constants  defined  to  be  the  three  terms  in  the 
numerator  of  B-54 


T1  -initial  value  of  spatial  coordinate  throughout  a 

given  mesh  interval 

T2-4  -intermediate  values  of  the  spatial  coordinate 
throughout  a  given  mesh  interval 

THETAE,  -Grande. 

THETAG 


TN 

X 

Y 

Z 


final  value  of  spatial  coordinate  where  numerical 
integration  ends 


normalized  electron  density. 


-normalized  negative  ion  density,  y*  A8 

w»* 


-final  estimate  of  w  (as  defined  in  3-6)  at  the  end  of 
each  iteration  in  the  integration  routine 


Z3 


-estimate  of  w  in  integration  routine 


1000  NU-14.1:  ALPHA* 7 

1010  MUE=1022.7:  MOI*2.25:  MUN-4.4 

1020  SE*( NO-ALPHA) /HOE:  SI«NO/MOI:  SN-ALPHA/MUN 

1030  THETAE-2.4:  THETAG-.15 

1040  C2* (SE+SI+SN)/ (THETAE+THETAG) 

1050  R*7 

1060  M»2*THETAG*R/ (THETAE+THETAG) :  NC*1 

1070  GAMMA-THETAE/THETAG 

1080  N*100 

1090  Tl*0 

1100  TN*1 

1110  Z3»M+NC 

1120  Z*Z3 

1130  Q1=0 

1140  Q*Q1 

1150  DT* ( TN-Tl ) /N 

1160  FOR  N*1  TO  100 

1170  Z3=Z 

1180  T4*T1+N*DT 

1190  T3*T4-DT/2 

1200  T2=T4-DT 

1210  K1=Q 

1220  GOSUB  1430 

1230  L1=-C2*X 

1240  Z3*Z+Kl*DT/2 

1250  K2*Q+Ll*DT/2 

1260  GOSUB  1430 

1270  L2»-C2*X 

1280  Z3*Z+K2*DT/2 

1290  K3*Q+L2*DT/2 

1300  GOSUB  1430 

1310  L3*-C2*X 

1320  Z3*Z+K3*DT 

1330  K4«Q+L3*DT 

1340  GOSUB  1430 

1350  L4»-C2*X 

1360  Z-Z+(Kl+2*K2+2*K3+K4) *DT/6 

1370  Q=Q+(Ll+2*L2+2*L3+L4)*DT/6 

1380  Z3=Z 
1390  GOSUB  1430 

1400  E=MUP*THETAG* ( 1 +GAMMA+2  * ALPH A*G AMMA) *Q/(MUE* (1+MUP 

* (1+ALPHA) /MUE+ALPHA*MUN/MUE) *X)  -  THETAE*Q/X 
1410  PRINT  "t*";T4,"e*";E,"X*";Xf "Y*" ;X* GAMMA 
1420  NEXT  N 

1430  REM  *SUBROUTINE— FIND  X~GAMMA* 

1440  LO*0 
1450  HI-1 
1460  X-.5 

1470  DIF«Z3-M*X*GAMMA-NC*X 

1480  IF  ABS(DIF) < .0001  THEN  GOTO  1530 


1490  IF  DIF<0  THEN  HI=X 
1500  IF  DIF>0  THEN  LO»X 
1510  X»(L0+HI)/2 
1520  GOTO  1470 
1530  RETURN 


The  symbols  appearing  in  the  following  code  have  the  same 
definitions  given  in  appendix  E  with  the  exception  of  the 
following  minor  changes  and  additions: 

AA,BB, 

CC,DD  -A,B,C  and  D  as  defined  in  B-51  and  B-53 
NN,MM,S  -n,m  and  s  as  defined  in  B-49  and  B-53 

u  replaces  x  and  s  replaces  f  in  the  definitions  of  HI,  LO, 
M  and  NC  where  the  last  two  are  now  given  by  3-2. 

The  flow  chart  is  also  similar  to  the  one  given  in  appendix 
E  where  the  changes  noted  above  must  be  made. 


10  THETAE=2.4:  THETAG-.15 
20  GAMMA=THETAE/THETAG 
30  MUE=1022.7:  MUI*2.25:  MUN-4.4 
40  R»10 
50  ND=14 

60  ALPHA=NU*MUN*R/ (MUE+MON*R) 

70  SE= (NU- ALPHA) /MUE :  SI*NU/MUI:  SN=ALPHA/MUN 

80  A=GAMMA* (1+1/GAMMA) *SE 

90  B*SI  +  SN  -  SE/GAMMA 

100  C=SI  -  SN  +SE 

110  D«2*SE*R 

120  F- (B-C- ( (B-C) *2  +  4*A*D)*.5)/(2*D) 

130  PRINT  F 

140  NN=C+D*F :  MM»B-D*F 
150  S-NN/MM 
160  PRINT  S 
170  BB»D/(MM-NN) 

180  AA»(1-BB*(1-F))/(1-F)aS 
190  CC«AA*F:  DD«1+BB*F 

200  PRINT  "A»"jAA,*B»";BBf "C»";CC,"D*";DD 
210  C2*(SE+SI+SN)/ (THETAE+THETAG) 

220  M»AA+2*R*CC/ (GAMMA+1) :  NC=BB+2*R*DD/ ( GAMMA+1 ) 
230  N-100 
240  T1»0 
250  TN-1 

260  Z3-M*(1-F)*S  +  NC*  (1**F) 

270  Z-Z3 
280  Ql*0 
290  Q*Q1 

300  DT“ (TN-T1) /N 

310  FOR  N«1  TO  100 

320  Z3«Z 

330  T4«T1+N*DT 

340  T3-T4-DT/2 

350  T2-T4-DT 

360  K1«Q 

370  GOSUB  570 

380  Ll»-C2* (AA*U*S  +  BB*U) 

390  Z3-Z+Kl*DT/2 

400  K2«Q+Ll*DT/2 

410  GOSUB  570 

420  L2— C2*(AA*U“S  +  BB*U) 

430  Z3«Z+K2*DT/2 

440  R3»Q+L2*DT/2 

450  GOSUB  570 

460  L3»-C2*(AA*UaS  +  BB*U) 

470  Z3-Z+K3*DT 

480  K4«Q+L3*DT 

490  GOSUB  570 

500  L4— C2*(AA*U*S  +  BB*U) 


510  Z»Z+(Kl+2*K2+2*K3+K4)*DT/6 

520  Q-Q+(Ll+2*L2+2*L3+L4)*DT/6 

530  Z3-Z 
540  GOSUB  570 

550  PRINT  "t-" ;T4, "X»" ;AA*0*S+BB*0/ "Y-" ;CC*0AS+DD*0, *0 

560  NEXT  N 

570  REM  *SUBROOTINE— FIND  0* 

580  LO-O 
590  HI-l-F 
600  0-.5 

610  DIF«Z3-M*0*S  -NC*C 

620  IF  ABS (DIF) <.0001  THEN  GOTO  6 70 

630  IF  DIF<0  THEN  HI-0 

640  IF  DIF>0  THEN  LOO 

650  0- (LO+HI) /2 

660  GOTO  610 

670  RETORN 
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